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ABSTRACT 

We  study  a  model  equation  for  the  elongation  of  filaments  or  sheets  of 
polymeric  liquids  under  the  influence  of  a  force  applied  to  the  ends. 
Mathematically  this  equation  has  the  form  of  a  nonlinear  Volterra  integro- 
differential  equation  with  the  kernel  given  by  a  finite  sum  of  exponentials. 
The  unknown  function  denotes  the  length  of  the  filament  or,  respectively,  the 
thickness  of  the  sheet.  We  study  the  equation  both  analytically  and 
numerically.  The  force  is  assumed  to  converge  to  zero  exponentially  as 
t  ->•  -oo  and  to  vanish  identically  after  a  f inice  time  tg.  It  is  shown  that 
under  this  condition  there  is  a  unique  solution  which  approaches  a  given  limit 
as  t  *  -oo;  moreover,  the  solution  also  has  a  limit  as  t  +  +°°.  a  numerical 
scheme  is  analyzed  and  convergence  uniformly  in  t  is  established. 

Particular  attention  is  paid  to  the  dependence  of  solutions  on  a  parameter 
U,  which  corresponds  to  a  Newtonian  contribution  to  the  viscosity.  It  is 
proved  that  solutions  converge  uniformly  in  t  as  u  +  0,  and  that  the 
convergence  of  the  numerical  scheme  is  also  uniform  in  p. 
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SIGNIFICANCE  AND  EXPLANATION 


The  evolution  of  the  shape  of  a  filament  or  a  sheet  of  a  polymeric  liquid 
subjected  to  an  external  force  f(t)  is  described  by  the  equation 

yy(t)  +  a(t-s)  ^  -  y(  s )  )ds  =  f(t)ya(t) 
y  \ s) 

where  y  denotes  the  length  of  the  filament  or  the  thickness  of  the  sheet, 
respectively,  and  y  is  a  Newtonian  contribution  to  the  viscosity  which  can 
either  be  positive  or  zero.  The  exponent  a  depends  on  the  physical 

situation  under  study. 


We  regard  the  length  at  t  =  -«  as  known  and  investigate  its  evolution.  It 
is  shown  that  for  a  physically  realistic  class  of  functions  f  there  exists  a 
unique  solution,  and  that  the  length  approaches  a  new  stationary  value  at 
t  =  oo  (which  is  in  general  different  (greater)  from  the  value  at  t  =  -«°). 
Numerical  calculations  are  performed  for  several  functions  f.  For  the  kernel 
a  we  choose  values  given  in  the  literature  for  polyethylene  at  150°C.  Our 

computations  show  in  particular  that  the  solutions  do  not  significantly  depend 

_2 

on  y  unless  y  exceed:,  10  000  Nm  sec. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  authors  of  this  report. 


A  NONLINEAR  VOLTERRA  INTEGRODIFFERENTIAL  EQUATION 
DESCRIBING  THE  STRETCHING  OF  POLYMERIC  LIQUIDS 


*  *  * 

P.  Markowich  and  M.  Renardy 

1 .  Introduction 

In  this  paper  we  consider  a  mathematical  model  describing  the  stretching 
of  a  filament  or  a  sheet  of  a  molten  polymer  under  a  prescribed  force  f. 
These  two  physical  situations  are  illustrated  by  the  following  diagrams: 


-f 

Our  model  is  based  on  the  following  physical  assumptions: 

(i)  The  polymer  satisfies  the  "rubberlike  liquid"  constitutive  relation 
[51. 

(ii)  The  strain  and  stress  tensors  are  independent  of  spatial  coordinates, 
and,  in  particular,  inertial  forces  are  neglected  (for  a  model  that 
includes  inertial  forces  see  [9]). 

(iii)  The  molten  polymer  is  incompressible. 

This  material  is  based  upon  work  supported  by  the  National  Science 
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Under  these  assumptions  the  problem  is  described  by  the  equation  (for  a 


derivation  see  [6],  [9]): 

(1.1)  py(t)  +  /t  a(t-s)f^— — —  y{  s )  Ids  =  f(t)ya(t),  -«  <  t  <  °° 

y  ( s) 

where  y  denotes  the  length  of  the  filament  or  the  thickness  of  the  sheet, 
resp.,  p  is  a  nonnegative  material  constant  modelling  a  Newtonian 
contribution  to  the  viscosity,  which  may  physically  come  from  a  solvent  or 
fractions  of  low  molecular  weight,  and  the  memory  kernel  a  has  the  form 

N  -Xu 

(1.2)  a(u)  =  l  K  e 

4=1  1 

with  positive  constants  K  ,  X  .  f  denotes  the  force  acting  on  the  ends  of 

X  X 

the  filament,  or  -f  denotes  the  force  acting  on  the  edges  of  the  sheet, 
resp.  The  exponent  a  is  2  for  the  filament  and  ^  for  the  sheet  (for  our 
mathematical  analysis,  we  assume  0  <  a  <  3).  The  difference  comes  from 
geometric  reasons:  If  f  were  denoting  the  force  per  unit  area,  a  would  be 
1  for  both  cases.  Due  to  the  incompressibility,  however,  the  area  on  which 
f  is  acting  depends  on  y. 

Although  this  has  no  significance  to  the  mathematical  analysis,  the 
physical  relevance  of  the  model  is  limited  to  f  >  0  for  the  filament  and 
f  <  0  for  the  sheet.  If  e.g.  one  attempts  to  compress  the  filament,  then 
buckling  rather  than  contraction  would  be  observed,  and  this  instability  is 
not  described  by  our  equation. 

A  problem  related  to  ours  was  investigated  by  Lodge,  McLeod  and  Nohel 
[6],  They  assume  y(t)  is  given  for  t  <  0,  it  is  nondecreasing  (which 
implies  but  does  not  follow  from  f  >0),  and  y(-°°)  =  1.  They  then  assume 


f  =  0  for  t  >  0  and  study  the  elastic  recovery.  For  a  class  of  kernels 
a  and  functions  F(y(t),  y(s))  under  the  integral,  which  include  those 
specified  above,  they  prove  the  existence  of  a  unique  solution  to  the  history 
value  problem,  which  is  nondecreasing  for  t  >  0  and  converges  to  a  limit 
y(<»)  >  1.  Their  proofs  rely  on  monotonicity  arguments,  and  they  also  prove 
that  the  solutions  depend  monotonically  on  the  prescribed  history  and  the 
parameter  y.  One  of  the  main  points  in  their  analysis  is  the  behavior  of 
solutions  near  y  =  0,  in  this  case  the  solutions  become  discontinuous  at 
t  =  0,  and  they  face  a  singular  perturbation  problem  with  a  boundary  layer. 

On  the  basis  of  these  results  Nevanlinna  [8]  used  an  implicit  first  order 
Euler-type  discretization  scheme  for  (1.1),  He  proved  that  this 
discretization  preserves  all  the  monotonicity  properties,  and  that  the  global 

V 

error  is  0(h  )  for  any  y  <  1  uniformly  with  respect  to  y  @  [0,y]  and 
t  e  [tQ,<»),  tQ  >  0.  It  was  not  shown,  however,  that  the  scheme  is  first  order 
accurate  uniformly  in  t  and  y,  i.e.  that  the  global  error  is  0(h). 

In  our  analysis,  we  prescribe  a  continuous  function  f(t),  which 

satisfies  lim  e  0tf(t)  =  0  for  some  c  >  0,  and  f  =  0  for  t  e  [t^,®). 

We  prove  that,  for  any  such  f,  problem  (1.1)  has  a  unique  solution  y(t) 

satisfying  lim  y(t)  =  1.  This  convergence  is  exponential,  moreover,  the 

solution  exists  globally  in  time,  and  converges  exponentially  to  a  constant 

y(<*>)  >0  as  t  ♦  ®,  more  precisely,  we  have  lim  e  crt  ( y ( t )  -  1)  = 

t  +—  oo 

lim  e<7t(y(t)  -  y(«))  =0.  This  holds  for  any  y  >  0.  The  solution  depends 
t  ♦+>*> 

oo 

continuously  on  y  in  a  norm  stronger  than  the  L  -norm  (more  specif ically , 


in  an  exponentially  weighted  L  -norm,  which  incorporates  the  asymptotic 
behavior  as  t  +  ±°°)  ,  even  at  y  =  0.  No  boundary  layer  occurs,  since  the 
solution  for  y  =  0  has  the  correct  asymptotic  behavior  as  t  +  ±«>.  Our 
proofs  are  mainly  based  on  the  implicit  function  theorem  and  Liapunov  function 
arguments . 

In  the  second  part  of  the  paper  we  discuss  the  computational  solution  of 
(1.1).  Like  Nevanlinna,  we  use  a  first  order  implicit  Euler-type 
discretization  with  uniform  mesh  size  h,  after  having  cut  the  interval 
[-•, 0]  at  t  .  In  the  convergence  proof,  we  use  a  discrete  analogue  of 

OO 

exponentially  weighted  l  -spaces  (infinite  sequences  converging 
exponentially  on  both  sides).  Choosing  a  space  with  an  exponential  weight 
given  by  e *  0  e >1 1 1 ,  0  <  c  <  0,  we  obtain  an  error  estimate  of  the  form  0(h) 

-clt-J 

+  o(e  )  in  the  norm  of  that  space,  moreover,  this  holds  uniformly  in 

y  e  [  0 , °°)  and  e  e  [0,eQ]  ,  cQ  <  a-  The  main  tool  in  the  proof  is  Keller's 
[3]  nonlinear  stability  concept. 

Our  numerical  results  imply  that  the  solution  y(t,u)  does  not  differ 
significantly  from  y(t,0)  on  [-00  ,  00]  if  y  is  smaller  than  a  certain 
fairly  large  number.  If  y  exceeds  this  number,  then  the  solutions  change 
considerably . 

The  paper  is  organized  as  follows:  In  chapter  2  we  present  the 
analytical  results,  chapter  3  concerns  the  discretization  procedure,  and  the 
computations  are  reported  in  chapter  4. 
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2.  Analysis  of  the  Continuous  Problem 
Solutions  for  Small  Forces 


N  -X.  u 

Let  us  consider  equation  (1.1),  where  0  <  a  <  3,  a(u)  *  \  K.e  1  , 

i=1  1 

and  n  >  0.  This  equation  can  be  reduced  to  a  system  of  ODE's  in  two  ways. 
We  set 


g.(t)  =  ft  K.e  i 
i  i 


-X.  (t-s> 


1 


y2(s) 


ds 


-X  (t-s) 

l . (t )  =  f  K  e  y( s ) ds 

1  *  — OO  1 


Then  (1.1)  reads 


y  =  - 1  (  I  f^y3  -  V  -  f(t)Ya) 


i=1 


(2.1) 


K. 

*i  =  -X.q.  +  -f 

y 


h.  =  -X.h.  +  Kv 
l  li  i 


If  we  set  y,  =  g.y2,  6.  =  — ,  we  obtain 
i  i  i  y 


y  =  ~  ~  l  l  (Yi  -  s.>  y  -  f(t)ya) 

p  i=  1 


(2.2) 


=  -XiY.  +  -  *  y.  I  (y.  -  6.)  +  \  y.f(t)y 


a-  1 


6.  =  -X.6.  +  K  +  —  8,  y  (y.  -  6.)  -  -  S.f(t)y 
l  li  i  U  i  1  1  pi 


a-1 


Both  forms  (2.1)  and  (2.2)  will  be  used  in  the  following. 


K. 

i 


Clearly,  if  f  =  0,  then  y  =  1 ,  qi  =  h^  =  —  =  y.  =  A  is  a  stationary 

i 

solution. 


-*>- 


LEMMA  2.1.  The  2N  +  1 -square  matrix  setting  up  the  right  hand  side  of  the 
linearization  of  (2.1)  (or  (2.2))  at  the  stationary  solution  y  =  1,  gt  =  Iv 


—  has  zero  as  a  simple  eigenvalue.  All  other  eigenvalues  have  negative 
i 

real  parts. 

PROOF :  Clearly,  (2.1)  and  (2.2)  give  the  same  eigenvalues.  Let  us  consider 


;  2. 1 ) .  The  linearization  is  set  up  by  the  following  matrix 


A  = 


N  3K 

-  y  — 


i=i  ^i 


2  2 
u  u 


-2K. 


-X. 


-2K„ 


0  “  X„ 


0  0  . 

0  0  . 


-2K 


0  0 


-X 


0  0 


.  0 


K. 


0  0 
0  0 


-X 


1 

0-  x„ 


.  0 
.  0 


0  .  .  . 


0  0 


.  .-X 


This  yields  the  characteristic  polynomial 


3K . 


3K . 


p ( x )  =  n  (-x.  -  x)  (-1  ~r  -  \  -  I  ,  , 

3  3  i  uXi  i  Li(  -  XA- X) 


Thus  N  eigenvalues  are  give  by  A  =  -A  ,  the  remaining  N  +  1  eigenvalues 


are  the  zeros  of  the  last  factor.  Obviously  one  of  these  is  zero,  and  it  is 
simple.  It  remains  to  be  proved  that  all  the  remaining  roots  have  negative 
real  parts.  Consider  the  equation 


(2.3) 


3K.  3K. 

_  y  _ i  _  x  _  v  1 

uX.  y(-X.-X) 

li  i  i 


=  0 


-6- 


and  its  sign  is  positive  for 


The  left  hand  side  has  poles  at  X  =  -X  , 

X  -*•  -X.+  and  negative  for  X  +  -  X.-.  For  convenience,  let  the  X. 's  be 
i  i  1 

ordered  such  that  X^  <  X^  <  •  •  •<  X^.  It  follows  that  there  is  a  root  in  each 

interval  (-X.,  -X.  .)  and  another  root  between  -X  and  Hence  all 

1  i+1  N 

non-zero  roots  are  real  and  negative. 

We  want  to  prove  the  existence  of  solutions  for  small  f  using  the 
implicit  function  theorem.  The  spaces  in  which  we  apply  this  theorem  are 
defined  in  the  following: 

DEFINITION  2.2.  Let  Y0,n  =  {g  6  Cn(R,  R)  |  lim  e°^g(k)(t)  =  0  for 

n  It  t  X” 

k  =  0,1,..., n}.  A  natural  norm  in  Y°'  is 


llgll  =  l  sup  |e°'t^gtk>(t)  | 
k=0  teR 

Moreover,  let  xa,n  =  {f  e  Cn(R,  R)  |  lim  e°^t^k'  (t)  *  0  for  k  =  1,...,n, 

|t  |  >« 


3  f(°°)  such  that  lim  e^tflt)  -  f(®>))  =  lim  e  ^ftt)  =  0} 
t  -H»  t  00 


Or  n 

A  natural  norm  in  X  is 


flf  II  =  J  sup  |  e  °  ^ ^  f  ^k  \  t )  |  +  sup  |e  crt  f  { t )  |  +  suple^ffft)  -  f  (  =>)  I 

k=i  teR  t<o  t>o 


+  I  f  ( °°)  I 

THEOREM  2.3.  Let  Y  denote  ( Y , I . f 1, < • • • r  Y„ »  5 , , 6, , . . . , 5„ )  and 
— -  .  12  N  I  2  N 

K1  kn  Ki  kn 

Yn  =  (l,  — Let  o  >  0  be  small  enough  (smaller  than  all 
X1  *N  A1 

the  absolute  values  of  the  non-zero  eigenvalues  of  A).  Then  the  following 
holds:  If  f  e  Y°'n  has  sufficient  small  norm,  then  (2.2)  has  a  solution 

Y  satisfying  Y  -  Yfl  e  X0,n+1  x  (Y0,n+1)2N.  Y  depends  smoothly  on  f. 


Proof :  When  we  put  Y  -  Y.  =  Z#  equation  (2.2)  can  be  written  in  the  form 


G(Z,f)  =  0,  and  G  is  a  smooth  mapping  from  (X 


x  (Y 


')*">  x  y0'1 


into  (Y°'n)^+\  Moreover,  the  linearization  DzG(0,0)  is  the  mapping 


N  2K . 

(y ,  Y,  ,  <5,  )  (~py  +  y  ( "Y .  —  6 .  > ,  y  +  \.y.  +  /  (  y  .  -  6 . ) ,  6.  +  X.5.  - 


i  i 


i=  1 


l  i 


yX.  f 
i  3 


y  (y  ~>5^>).  According  to  lemma  2.2,  the  y  and  6  components  form  an 
isomorphism  from  (y0,n+1)2N  onto  (Y°'n)2N.  Moreover,  the  mapping  y  y 
is  a  bijection  from  Xa,n+1  on  Y0,n.  Therefore  DzG(0,0)  is  an  isomorphism 
from  X°'n+1  x  (ya,n)2N  onto  (Ya'n)2N  +  \  The  implicit  function  theorem 


yields  the  result. 
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Global  Behavior  of  Solutions  for  Large  f 


Theorem  2.4. 

Let  u  >  0 

and 

f:  R  +  R  be  continuous  and  such 

that 

lim 

t  +-  “> 

e  ^fft) 

II 

O 

o 

V 

o 

is  as 

in  Theorem  2.3),  f(t)  =  0  for 

t  >  t  .  For 

0 

every 

such  f , 

equation 

(2.2) 

has  a  unique  solution  satisfying 

lim  y(t)  =  1 
t+- 00 

K. 

lim 

t+-°° 

y.  =  lim 

1  t-*~ 

<5.  =  r1. 
l  X. 

00  1 

This 

solution  exists  globally  for  all 

times  t. 

and 

lim  y(t) 

exists  and  is 

strictly  positive. 

t 

Proof :  If  is  chosen  large  enough,  e  becomes  small  on  (-°°,-t  ) 

and  one  can  use  an  implicit  function  argument  analogous  to  theorem  2.3  to 
prove  the  existence  of  a  solution  on  (-«,-t  ).  This  solution  is  unique  in 
the  class  of  solutions  approaching  their  limiting  values  at  t  =  -»  at  a  rate 
of  e0^.  However,  if  a  solution  tends  to  these  limits  at  all,  it  can  be  seen 
from  the  last  two  equations  of  (2.2)  and  the  implicit  function  theorem  that 
Y .  and  6^  tend  to  their  limiting  values  at  a  rate  of  e0^.  The  first 
equation  then  implies  that  y  approaches  its  limiting  value  at  the  same 
rate.  Hence  the  solution  is  actually  unique  in  the  class  of  all  solutions 
approaching  the  prescribed  limits  at  t  *  as  claimed  in  the  theorem. 

We  now  continue  this  solution  to  the  right,  and  we  have  to  make  sure  that 
it  does  not  blow  up  at  a  finite  time.  For  that  purpose  it  is  more  convenient 
to  consider  (2.1)  rather  than  (2.2).  From  the  second  and  third  equation  we 
see  that  as  long  as  y  stays  positive,  q^  and  h^  have  a  positive  lower 

bound  for  all  finite  times,  which  is  independent  of  y.  Hence,  if  y  becomes 
3  a 

too  larae,  q^y  will  dominate  over  fy  and  also  over  h^  (since  this  is 

less  than  some  constant  times  max  y(x)).  Analogously,  if  y  becomes  too 

(  -°°,t] 


smalL,  h.  will  dominate  over  fya  and  g^y^.  Hence  v  cannot  go  to  zero 
or  infinity  in  finite  time,  whence  we  find  global  existence. 


_g_ 


Let  now  t  >  tp.  Then  f  =  0,  and  using  (2.2)  again,  we  find  the 
Liapunov  function 


(2.4) 


N 

V 

L 

i=1 


a.  a.  3.  8. 

11  i  i  ii 

2  K.  M  K. 


a.  +  — 

i  X, 


si +  r 


N 

~  I 

i=  1 


X.  u  a, 

l  l 


K~  +  V 


a.  +  “r~ 
i  X. 
i 


ei  +  r 


l  (a.  -  8.  ) 

''ii 

i=1 


K.  K. 

Here  we  have  put  a,  =  Y.  -  -r~»  6,  =  6.  -  — •  As  we  know  that  v.  and  6. 

i  iX  iiX  l  l 

i  l 

K.  K. 

stay  positive,  the  denominators  a.  +  — ,  8.  +  t—  are  always  positive,  and 

1  A,  1  A . 

1  1 

the  left  side  of  the  equation  (2.4)  is  thus  the  derivative  of  a  positive 
definite  function  that  decreases  along  trajectories.  As  an  immediate 
consequence  we  obtain  that  <x  and  8^  tend  to  0  exponentially  for 
t  -*■  <*>.  One  easily  concludes  from  (2.2)  that  In  y  approaches  a  constant. 


and  hence  lim  y(t)  >0  exists. 

t-K» 


The  next  corollary  provides  information  on  the  final  recovery  for 
physically  significant  forces. 

Corollary  2.5.  If  f  is  always  non-negative  and  not  identically  zero,  then 
y(°»)  >  y(-°°),  if  f  is  always  non-positive  and  not  identically  zero,  then 
y(«>)  <  y(-“). 

Proof :  Assume  f  >  0,  the  other  case  is  analogous.  It  is  immediate  from  the 

integral  equation  (1.1)  that  f  >  0  implies  y  >  1  for  all  t.  Moreover,  if 

*  ★ 

f  /  0,  there  must  be  some  t  such  that  y(t  )  >  1.  Let  now  z(t)  = 


mm 


y(T).  Then  (1.1)  implies  that 


z(t)  >  min  (0,-  -  /fc  a(t-s) (z3(t)-1)ds) 

dt  +  y  J-oo 

* 

>  -  —  a<  t-s)  <  z3(  t )  —  1  )ds 

If  z(t)  -  1  is  sufficiently  small,  this  gives  an  inequality  of  the  form 

[j-l  z(t)  >  -  Ce  kt(z-1)  . 

at  + 

It  follows  immediately  that  lim  z(t)  >  1. 

Remark:  These  results  are  obviously  expected  on  a  physical  basis.  Namely, 
they  simply  state  that  pulling  the  filament  effectively  increases  its  length 
(f  >  0,  a  =  2)  or  the  thickness  of  the  sheet  decreases  (f  <0,  a  =  ^) , 
resp. 

We  now  give  an  argument  showing  that  theorem  2.4  does  not  hold,  if  the 

condition  that  f(t)  =  0  for  t  >  tQ  is  replaced  by  exponential  decrease 

of  f  and  a  i-  1  (in  case  a  =  1  the  previous  argument  still  goes  through, 

N 

the  only  difference  being  that  f(t)  \  (a.  -  8. )  has  to  be  added  on  the 

i=1 

right  side  of  (3.1)).  We  restrict  ourselves  to  the  case  N  =  1.  (2.1)  reads 

•  3  ci 

-uy  =  gy  -  h  -  f(t)y 

•  ,  K 
g  =  -Xg  +  — 

y 

h  =  -Xh  +  Ky 

We  solve  these  equations  for  t  >  0  by  the  following  ansatz: 


vt 

-2vt  -Xt  , 

* 

1 

ft 

y  =  y0e  >  <?  = 

v  +  gie  '  h  = 

nQe  +  ^e 

^  (l-a)vt 

3-a  ( (  3-  a)  v-  X)t 

.  -a  (-av-X)t 

f  =  foe 

+  Vo  e 

-  Vo  e 
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After  some  calculation  one  finds  that  this  satisfies  the  equations  if 


30  2 

Y0  ( X-2 v) 


Kyo 

'  ho  =  and 


foyo 


a-1  _  3vK  t  \iv(  1-2  v)  (  It  y) 
(  X—  2  v )  (  X+v) 


We  thus  find  solutions  where  f  qoes  to  zero  exponentially,  but  y  ♦  °>  fcr 
a  >  1  and  y  -*•  0  for  a  <  1 . 

All  we  have  to  make  sure  is  that  by  appropriate  continuation  for  t  <  0 
we  can  match  the  conditions  at  t  =  -<*>.  For  this  purpose  continue  y  in  an 
arbitrary  way  to  the  left  such  that  y  is  smooth  and  approaches  1 
exponentially  at  t  =  The  equations  for  g  and  h  then  have  unique 

solutions  approaching  —  for  t  *  These  solutions  can  be  matched  to  the 

A 

solutions  for  t  >  0  by  appropriate  choice  of  g1  and  h^  Finally  f  is 
determined  by  the  first  equation. 
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The  Case 


In  this  case  the  first  equation  of  (2.1)  becomes 


y3  *  I  9.  -  l  h  -  f(t)ya  =  0 
i=1  i=1 


Proposition  2.6.  For  any  g  >  0,  h  >  0  and  0  <  a  <  3  the  equation  F(y) 
gv3  -  h  -  fya  =  0  has  a  unique  solution  in  (0,“). 

Proof :  We  have  F(0)  <  0,  lim  F(y)  >  0,  so  there  is  clearly  a  positive 

y  4co 

solution.  To  show  it  is  unique,  we  investigate  zeros  of  F' (y) .  We  have 

F'(y)  =  3gy^  -  ofya  \  If  y  >  0  and  F'(y)  =  0,  we  find  F(y)  = 

1  3  3 

—  yF'(y)  +  y  (1  -  — )g  -  h  <  0.  This  means  F  cannot  have  a  positive 
a  a 

maximum,  whence  the  result. 

The  solution  y(g,h,f)  can  then  be  inserted  into  the  other  equations, 
yielding  a  system  of  2N  equations. 

Theorem  2.7.  The  same  statement  as  in  Theorem  2.4  holds  also  for  y  =  0. 
Also,  Corollary  2.5  still  holds. 

Sketch  of  the  Proof :  The  existence  of  a  solution  on  (-“,-t  )  and  global 
existence  in  time  are  proved  in  the  same  manner  as  before,  and  we  do  not 
repeat  the  arguments.  If  f  =  0,  one  finds  from  (2.2) 


y.  =  -X.  y.  +  K.  +  2y.  x 
l  li  l  i  y 


6.  =  -X,  6.  +  K.  -  5.  L 
i  li  i  l  y 


This  leads  to 


N  a.  a.  B.  (!.  N  X.  a.  X.  B.  •  N 

r  .  1  11  11,  r  11  11,  y  r 

.'■  2  K,  K,  1  '■  K.  K,  y  i 

i=1  1  1  1=1  „ ,  1 ,  1  i=1 

a.  +  t~  B.  +  7—  2(a,  +  — )  B.  +  T~ 

1  X.  1  X.  1  X.  1  X. 

11  11 


where  a,  and  B.  are  defined  as  before. 
1  1 
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r 


Since  £(ol  -  8^)  is  now  equal  to  zero,  we  still  find  that  and  0^ 

approach  0  exponentially,  whence  the  result. 

For  the  corollary,  observe  that 


y(t)  »3y2(t)  a(t-s)  a1  (t-s)  [—~  —  -  y(s)]ds  . 

y  (s)  y  (s) 

Using  this,  one  can  apply  an  ergument  analogous  to  the  previous  one. 

Finally,  we  want  to  prove  that  solutions  depend  continuously  on  y,  even 
at  vi  =  0.  Let  f  e  ya'n  be  given  such  that  it  either  has  a  small  norm  or  it 
satisfies  the  conditions  of  theorem  2.4.  We  know  that  a  unique  solution  y(t) 
satisfying  y(-°°)  =  1  exists  both  for  y  =  0  and  for  u  >  0.  In  (2.1),  we 


N  N 

put  g  =  l  g^,  h  =  [  h^,  and  z  =  y 


i=1 


i=1 


o  /Vi 

—  (for  p  =  0,  f  =  0,  the  first 


equation  of  (2.1)  is  solved  by  y  =  We  obtain 

•  (-3  /h  \3  .  -,fcv/-3/h^  >,0  .  d  3  /h 

-uz  =  g(  yf- *  z)  -  h  -  f(t>(y-  +  z)  +  ii  ^  y- 


gi  = 


"  Xigi 


(2.5) 


+  TTTh  \2 

(vr  z> 


h.  - 

l 


A.h  +  K.  (t/~  +  z  ) 
i  i  V  g  ‘ 


As  we  have  proved,  there  exists  some  y^  >  0  such  that  for  every  y  6  [ 0 ,  vj^  1 


system  (2.5)  has  a  unique  solution  in  the  Banach  manifold 


M  =  { (z,  g  ,h  )  |  z  e  Y 
n  li 


o,n 


gi  “ 


\  3  /S_ 
i  V  2 


1C  — 

„  0,n  i  3  /h  „a,n, 

Y  '  hi  -  “  Vg  e  Y  }-  In 


V  g 

particular,  let  zQ,  gi  g,  h^  0  denote  the  solution  for  y  =  0. 

Linearizing  at  this  solution  (or  likewise  at  any  solution  for  y  >  0), 
we  obtain  a  system  of  linear  ODE’s  with  a  matrix  approaching  a  constant  limit 
as  t  ♦  and  t  +®.  From  a  discussion  of  the  asymptotic  behavior  of 


solutions  of  the  linearized  system  for  t  ±®, 


one  can  easily  see  that  for 


any  inhomogeneity  in  (Ya'n)/N+1  there  is  a  unique  solution  in  the  tangent 
space  of  Mn«  The  argument  parallels  our  existence  proof  for  solutions: 

First  consider  the  problem  on  f  — T)  with  T  large,  where  the  matrix  is 
approximated  by  the  linearization  at  the  trivial  solution.  Continuation  of 
solutions  for  t  >  -T  presents  no  problem,  since  the  equation  is  linear,  and 
finally  the  behavior  for  t  +  +»  must  be  discussed.  We  leave  the  details  of 
the  analysis  to  the  reader. 

Thus  the  linearization  is  a  densely  defined  bijective  operator  from  the 
tangent  space  of  Mn  onto  (y°'n)^N+\  It  is  thus  natural  to  attempt  proving 
the  existence  of  a  continuous  family  of  solutions  in  a  neighborhood  of  y  =  0 
using  the  implicit  function  theorem.  One  does,  however,  face  the  problem  that 
the  term  yz  represents  an  unbounded  operator. 

The  first  equation  of  (2.4)  has  the  form 

-y  3—  (z-z„)  =  p(t)(z-z„)  -  f (t) *L(h-h„ ,f-f )  +  nonlinear  terms  +0(y) 
dt  0  0  0  0 


where 


3  /h  3  fh 

p(t)  »  390  (  V  7"  +  2o)  "  af(t>(  V  + 


is  positive  and 


N  K 

converges  to  }  —  for  t  ±®.  L  is  linear  in  its  arguments,  and  the  term 

i=  1  1 

0  (  y)  does  not  involve  any  unbounded  operators,  after  the  second  and  third 


•  • 

equation  of  (2.4)  have  been  substituted  into  the  first  to  replace  g  and  h. 

It  is  easy  to  show  that  the  operator  ( y  +  p(t))  1  :  y°'n  -*■  y°'n  is 

dt 

strongly  continuous  with  respect  to  y.  Denoting  V  =  (z-z0,g1~g1  q,..., 
gN~gN  Q,h1-h1  g,...,hN-hN  n>,  we  can  thus  rewrite  (2.5)  in  the  abstract  form. 


(2.6)  L(y)V  =  N ( y , V )  <==>  V  -  (L( y) ) ~ ’n( y,V)  =  0 

where  L(y)  has  a  strongly  continuous  inverse  and  N(0,0)  =  0,  DvN(0,0)  =  0. 

The  existence  of  a  continuous  solution  V(y)  now  follows  from  the 
following  theorem. 
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Theorem  2.8; 


Let  X,Y  and  Z  be  Banach  spaces,  U  a  neighborhood  of  (0,0)  in 
XxY,  and  F;  u+Z  a  mapping  having  the  following  properties: 

(i)  F( 0, 0 )  =  0 

(ii)  F  is  continuous 

(iii)  F  is  continuously  differentiable  with  respect  to  y  for  each  fixed 

x. 

(iv)  DyF(0,0):  Y+Z  is  an  isomorphism. 

(v)  DyF  is  continuous  at  the  point  (0,0). 

Then  the  equation  F(x,y)  =  0  has  a  unique  resolution  y  =  f(x)  in  some 

neighborhood  of  (0,0),  and  f  is  continuous. 

The  proof  of  this  theorem  differs  by  no  means  from  the  standard  proof  of 

the  implicit  function  theorem  (cf.  [10],  [11]),  but  it  is  crucial  for  our 

problem  that  (iii)  and  (v)  are  sufficient  rather  than  continuity  of  DyF  in  a 

neighborhood  of  (0,0)  as  usually  required.  Namely,  we  can  identify  X 

with  R,  Y  with  the  tangent  space  of  Mn»  2  with  Y0'n,  x  with  y  and  y 

with  V.  For  y  fixed,  the  term  L( y)  ^N(y,V)  depends  smoothly  on  V, 

moreover,  since  lim  D  N(y,V)  =  0,  we  also  have 
y+0,V+0 

lim  Dy(L(y)  *N(y,V))  =  lim  L( y)  ^DyN(y,V)  =  0.  Hence  Theorem  2.8 
y+h,V+0  y-K),V+0 

applies  to  (2.6),  although  the  standard  form  of  the  implicit  function  theorem 

would  not.  This  yields  a  continuous  solution  V  =  V(y). 

d  -1  k 

Moreover,  the  mapping  (y,z)  «■  (y  —  +  p(t))  z  is  a  C  -mapping  from 
R  x  y°'n  into  ya,n  From  the  following  theorem,  which  was  also  proved  in 

[10],  [11],  one  concludes  that  V(y)  is  actually  a  C  -function  of  y  when 

.  .  ,  .  .  a,n-k 

regarded  as  lying  in  Y 
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Thgorem  2.9: 

Let  and  z^  resp.  (k=0,1,...N)  be  two  hierarchies  of  Banach 

spaces  such  that  Y^  c  Y^k+1^,  Z^  c  z^k+1^,  the  imbeddings  being 
continuous.  Let  X  be  a  finite  dimensional  Banach  space  and  P  a  mapping 
from  a  neighborhood  U  of  0  in  XXY^  into  Z*N^  having  the  following 
properties : 

(i)  F(U  n  (XxY(k)))  c  Z(k>  k=0 ,  1 , . . . ,N 

(ii)  For  each  fixed  k,  F,  :=F|  ,,  v  satisfies  the  conditions  of 

k  Un(X*Y(!c)> 

( k) 

Theorem  2.8,  when  it  is  considered  as  a  mapping  from  XkY  into 

l  v) 

Z'  .  For  x  fixed,  F^(x,.)  is  a  smooth  (i.e.  sufficiently  often 
differentiable)  mapping. 

(iii)  F:  XxY*k^  +  z(k+m)  ig  c^ass  for  each  k=0,1,...N  and 

m  <  N-k. 

(iv)  The  mapping  ( x,y  .u1 , . . .  ,u^  )  +  z  <=  D  i  j  F( x,y) (u1 , . . . ,u^ )  is 

K  y 

continuous  from  X  x  Y^  x  (Y^)"*  into  £  ^(X,Y^k+^ ) . 

The.1  the  solution  y  =  f(x)  e  Y^®  ^  existing  by  theorem  2.8  is  a  Cm-function 
of  x  in  some  neighborhood  Vm  of  0,  if  y  is  regarded  as  an  element  of 
Y(m>. 

We  summarize  our  results  in  the  following: 

Theorem  2 . 10: 

Let  f  e  Y0,n  be  given  such  that  either  f  has  small  norm  or  f(t)  r_  0 

for  t  greater  than  some  t  <  °°.  Then,  for  each  ye  [0,«>1,  (1.1)  has  a 

unique  solution  y  satisfying  y  -  1  e  X0,n.  In  the  limit  u  +  0, 

y  -  1  e  Xa,n  depends  continuously  on  u,  and  it  is  a  C'-funrtion  of  n 

when  regarded  as  dwelling  in  X0'0  k. 
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3.  The  Discretization  Scheme 

When  solving  (1.1)  numerically,  one  faces  the  problem  that  it  is  to  be 
solved  on  an  infinite  interval.  A  reasonable  way  of  doing  this  is  to  cut  at 


-T  <<  0,  and  replace  y(t)  for  t  <  -T 
thus  obtain  the  approximating  problem 


uy_T  +  /_*  a( t-s) ds  •  (y^T(t)-1) 


(3.1) 


by  its  limit  lim  y(t)  =  1.  We 
t+-<*> 


rt  ,y-T(t) 

+  J_T  a(t-s)  (— - y_T<.s))ds 

y  m(s) 

-T 


-  f(t)y“T(t)  =  0 


(3.2) 


y_T(t)  -  1  =  0,  t  <  -T  . 


On  the  finite  interval  the  integrodif ferential  equation  can  now  be  discretized 
in  a  straightforward  manner.  Like  Nevanlinna  [8] ,  we  use  a  first  order 
implicit  (Euler-type)  method,  because  for  this  simple  procedure  we  can  prove 
that  the  qualitative  properties  of  solutions  of  (1.1),  such  as  exponential 
decay  at  infinity  and  uniform  convergence  as  y  +  0 ,  carry  over  to  the 
discrete  problem.  Since  these  properties  are  essential  for  the  continuous 
problems  it  is  very  reasonable  to  require  that  the  computed  approximating 
solutions  exhibit  them  too.  Our  computations  have  shown  that  good 
approximations  can  be  obtained  with  quite  large  mesh  sizes,  and  so  the 
computational  effort  for  the  first  order  scheme  remains  reasonably  small. 


We  choose  me3h  points  t^  =  ih,  i  e  Z,  where  t_m  =  -T,  and  denote  by 
the  approximation  to  y(t^)  (or,  respectively,  to  y_T(t^)).  Then  our 
discretized  form  of  the  equation  reads 

3 

y  -y  t  i  y . 

U  —  a(t.-s)ds  •  (Yj-1 )  +  h  l  a  ( (  i—  j  )  h  )  (~  -  y  )  - 

j=-m+1  yT  3 

(3.3/  J  3 


-  f  ( t .  )  y  ?  =  0  ,  i>-m 

1  1 
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(3.4) 


1 


0 


i  <  -m 


y 


i 


Obviously, 


(3.5) 


/  m  a(t.-s)ds 

j  —  oo  i 


N 

l  t* 

i=  i  i 


V^nfV 


Equation  (3.3)  has  the  form 

(3.6)  Ciyi  +  ^  +  Vi  =  C4  ' 

where  the  c’s  depend  on  u,  h,  t_ffl,  t^  and  y.. ,  j  <  i.* 

The  analysis  of  the  discrete  equation  will  be  carried  out  in  the  same 

sort  of  spaces  as  the  analysis  of  the  continuous  equation.  We  therefore 
define  discrete  analogues  of  the  exponentially  weighted  spaces  introduced  in 
Definition  2.2. 

Definition  3.1: 

Let  X°:  =  {f  «  ({  )m  e  £°°|  lim  f  =:  f  exists, 

h  i  i=-<*>  l  °° 

i-H» 

lim  e^a^|f  -f  |  =  0,  lim  e  X°^|f.|  =  0} 
i>«  1  *  i+-„  1 


and  Y^:  =  {g=(g.).  e  SI  |  lim  e*  *  |  g .  |  =  lim  e  ioh  |  g,  |  =  0}  .  The 

n  i  i=— oo  i  i 

1  -foo  1  00 

natural  norms  in  these  spaces  are 

H?  II  =  sup  e^°^|f  . -f  |  +  sup  e  ^  ^  |  f  |  +  |f  | 

O  1  oo  1  OO 

x°  i>o  icn 

n 


and 


Setting 


2 


(y,-i>. 

1  1=-00 


l  oh  i  )  ”  i-  oh  |  i 

sup  e  | g . I  +  sup  e  | g . | 

i>0  1  i<0  1 

we  rewrite  (3.3)  and  (3.4)  in  the  abstract  form 


Equation  (3.6)  will  be  discussed  at  the  end  of  this  chapter. 
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(3.7) 


Ph,m<z>  =  0 


Now  let  f  e  Y°'n(o  >  0,  n  6  N)  be  given  such  that  the  assumptions  of  theorem 
2.4  hold.  It  is  an  easy  exercise  to  show  that,  for  any  e  e  [0,a), 

(3.8)  F. 

(we  explain  below  why  e  is  introduced). 


:  X0~£  -  Y°~e 
h,m  h  h 


The  aim  of  the  following  analysis  is  to  prove  that  (y^)^_  ^  converges 

00  Ct~  £. 

to  (y(t.)).  in  the  topology  of  X,  .  The  proof  will  be  based  on 

l  i=-»  h 

Keller's  [3]  nonlinear  stability-consistency  concept. 

Let  us  first  show  consistency.  The  local  discretization  error 
00 

£,=  (£.).  is  defined  by 

l  i=-°° 


(3.9) 


Z  =  F  ( (y(t  )  -  1)  ) 

h,m  1  i=-oo 


For  f  e  Y0'1  (which  implies  \sy  e  Y0'2,  uniformly  in  p),  we  find,  using 

the  exponential  decay  of  y  as  t  ->■  +«>,  that  in  the  limit  t  *  h  0 

-m 


-0|t  | 

(3.10)  (a)  U|  <  o( 1 )e  ,  i  <  -m  , 


Y(t  )  -  y(t  )  -alt  | 

(b)  In  - - - yy’  (t.  )  |  =  const.  o(  1  )he 

h  i 


( c) 


t_  t  y3( t . ) 

f  m  a(t.-s)ds  •  (y  (t.)  -  1)  -  f  a(t.-s)[— — — 

J  —  oo  1  1  J  —  00  iv2 

y  ( s) 


-  y ( s ) ) ds | 


<  const.  o( 1 )e 


o(2t  -t.) 
-m  l 


(d) 


i  y(t.)  t.  y  ( t . ) 

I  h  J  a(  (  i- j  )h)  (— - -  y(  t  . )  )  -  /  1  a(t.-s)(— - —  -y(s))ds 

j=-m+1  y  (t.)  -m  1  v  (s) 


-alt. 


<  const.  o{ 1 )he 


Here  o(  1  )  stands  for  a  factor  that  vanishes  as  t.  ♦  -°°.  Therefore, 

l 
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-  ol  t  |  -ol  t  -2t  I 
,  1  1  “in 

he  +  e  ,  i  >  -m 


(3.11)  |  (  <  const.  o(1) 


■al  t , 


,  i  <  -m 


From  definition  3.1  we  conclude  that 


-e|t  I 

(3.12)  Ul  <  const,  (h  +  o(e  "  ))  . 

ce 

The  constant  is  independent  of  h,  t_m,  0  <  y  r  °°,  0  <  e  <  eQ  <  a.  Note  that, 

- e I t_m | 

in  particular,  the  error  estimate  contains  a  term  o(e  ).  The  reason 

for  this  is  that,  when  approximating  (1.1)  by  (3.1),  (3.2),  we  have  replaced 

f  by  0  for  t  <  -T,  and  in  the  norm  of  ycre'n  this  introduces  an  error 
-e|T| 

of  the  order  o(e  ).  This  is  the  reason  why  we  have  introduced  the  e; 
for  e  =  0  we  would  still  get  convergence,  but  no  estimate  for  the  order. 

(3.12)  settles  consistency. 

For  the  stability  analysis,  we  calculate  the  Frechet  derivative  of  F^  m 

00 

at  the  exact  solution  (y(t.)  -  1)._  ,  which  is  denoted  by 


(3.13) 


L.  :  =  D„F  ( (  y(  t .  )  -  1)  .  ) 

h,m  z  h,m  i  i=-°» 


For  u  =  (u  ) .  @  X,  we  obtain 

i  i=-oo  h 
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<3*14' 


u^,  i  <  -m 


u  -u.  ,  t 

1  l- 1  .  r  -m 


3  J  a( t . -s) ds  y( t  )  u 

j-oo  i  i  i 


i  y2(t  > 

+  3h(  £  a(i-j)h)  — - )u  . 

j=-m+1  y  ( t . )  1 

3 


i  y  <t  > 

-  h  l  a(  (  i-j )h)  (2  — — —  +  1)u. 


j=-m+1 


y3(t.) 


orl 


-  a  f (t. )  y(t. )  u,  ,  i  >  -m 

i  i  1 


Stability  means  that  exists  and  that  it  is  bounded  as  an  operator  from 

n  ,n\ 

<  e  <  cQ  <  o. 

e  ?r*. 

h 


2, 

y  ( t .  i 


yfe  into  x°“e 
h  h 

uniformly  with 

respect  to 

h'  W 

U  and 

Therefore 

we  look 

at  the  equation 

A  A 

Vm0  =  v 

for  v  = 

00 

(v.  )  . 
l  i= 

For 

i  <  -m 

we  find  u^  = 

,  and  for 

i  >  -m 

we  show 

t 

G  (h,-m,ia)  =  ^  +  3  J 

a(ti~s)y2(t 

, ) ds  +  3h 

l 

i-1 

I  ' 

(3.15) 

a. 

ii  - 

i=-m+ 1 

- 

a  f  ( t .  )  v 
i 

a"  1  ( t .  ) 
l 

v  (t . ) 
3 


which  is  the  coefficient  of  in  (3.14),  is  bounded  away  from  0 

uniformly  in  h,  t  and  u. 

It  is  easy  to  show  that 


(3.16) 


G  ^  (  h ,  -m ,  y) 


=  ^  + 
h 


2, 

y  (ti) 

2 

y  ( s) 


ds 


a  f  ( t .  )  v 
i 


a- 1 


(t.  ) 
l 


-o(  It.  Mt  I) 

+  0(h)  +  o(e  1  "IT’  )  . 
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Using  (1.1)  we  get 


G.(h,-m,U)  =  £  + 


t . 


(f  1  a( t . -s)y( s ) ds  -  uy(t.) 


(3.17) 


i  i 

-o(  1 1.  |  +  |  t  | ) 


+  (1  -  y)f(t)ya(t.)  )  +  0(h)  +  o(e  1  _m  )  . 


It  follows  from  chapter  2  that 


(3.18) 


o  <  y0  <  y(t,p)  <  yQ,  !y(t,u)|  <  y 


uniformly  for  u  e  For  f  <  0,  (3.16)  provides  a  uniform  lower  bound 

for  G^ ,  and,  for  f  >  0,  (3.18)  provides  a  uniform  lower  bound,  since 

a  <  3. 

The  preceding  considerations  make  it  apparent,  why  the  term 
t_m  3 

f  a(t.-s)ds(y.  -  1)  should  be  maintained  in  (3.3).  If  this  term  were 

neglected,  the  uniform  bounds  on  would  no  longer  hold,  and,  unless  a 

constraint  of  the  form  ^  >  const,  is  imposed,  an  artificial  boundary  layer 

h 

can  be  venerated  at  t  We  see  from  the  above  that  L,  _  can  formally  be 

•  -m  n,m 

inverted.  It  remains  to  be  proved  that  the  solution  u  of 
(3.18) 

satisfies  an  estimate 

(3.20)  Hull  <  const.  Ilv  II 

,0'E  Y0_e 

\  h 


Lh,mu 


A  A 

V 


with  the  constant  independent  of  h,  t  ,  u  and  e. 

We  begin  with  the  reduced  problem  for  u  =0.  (3.14)  yields 


( 3.21 ) 


ir1  0 

u.  =  h  )  a ( ( i- j ) h ) a .  .  (h,-m)u.  +  v. 

i  .  '  ,  i -3  1  i 

l=-m+  1 
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where 


(3.22) 


a.  .  (h,-m) 
0  3 


y3(t  ) 

2  -3 -  +  i 

y  (t.) 


G.(h,-m,0)  '  V1  G.(h,-m,0) 

l  l 


with  Gi  as  defined  in  (3.17).  Since  y{-<=°)  =  1,  this  implies,  after  a 
simple  calculation 


a.  . (h,-m) 
i.3 


N  K. 


+  y  =  :  C(h, y) 


h  I 


.  .  X  „h 
11=1  l  . 

e  - 1 


where  y  -*■  0  as  h  -*■  0,  t  ♦  -<»,  t,  +  -“. 

-m  l 

Using  the  form  of  a,  we  get  from  (3.21) 


i-1  N  -X  (t.-t.) 

(3.24)  |u. |  <  C(h,y)h  £  £  K  e  1  X  1 

1  j=-m+ 1  Z=1  * 


u  I  +  Iv.| 


The  solution  w.^  of  the  equation  obtained  by  replacing  |v^|  by  the 
i( o-e )h 


larger  quantity  e 


llvll  ,  where 
00 
l 


(3.25) 


i  (  a-  e )  h  -  -  -  «  _  . o 

v .  =  e  v.  ,  v  =  (v.  )  .  el 

l  ll  i=-°° 


provides  an  upper  bound  for  I u^ | .  In  analogy  to  chapter  2,  we  substitute 


( 3.25)  g*  =  K^h  Y 

j=-m+1 


i-1  -X  (t.-t.) 

e  Z  1  1  w.  ,  l  =  1 ,2 . N  , 

1 


which  yields  the  difference  equation 
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(3.27) 


»i+1  -  gi  "  V  Wi  +  (e  ■’>«{ 


and  the  following  difference  equation  for 


N  -X  h  N  -X  h 

wi+1  -  wL  =  C(h,y)  l  (e  -1 )g^  +  C(h,Y)  l  K£he  1  wi 
£=1  1=1 


(3.28) 


.  r  (o-e)h  i  i(cr-e)h  - 
+  (e  -1  j  e  Dv  II 


(3.27)  and  (3.28)  form  a  system  of  difference  equations.  Setting 


(3.29)  z  =  (w  ,g\...,g?>#  <p  *  ((e(0_G)h  -  1)  ei(o_£)h  Ivl  ,0, 
l  i  i  i  i  « 

£ 


,0) 


we  can  rewrite  this  system  in  the  form 

(3.30) 
where 


zi+1  =  (I  +  A(h,y))zi  +  ^  i  >  -m 


N  -X  h  -X  h 

C(h,Y)  l  K  he  1  C(h,Y)  (e  -  1) - C(h,Y)  (e 


V 


£=  1 


(3.31)  A(h, y>  = 


K^he 


'X,h 


■V-  1 


Ve 


■v 


•  -V*  1 

e  -  1 
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The  solution  of  (3.30)  is  given  by 


(3.32)  z  =  (I  +  A(h,Y))1~1+mz  .  ,  +  l  <1  +  A<h,y))1  3 

i  -m-r  i  ,  j 

with  the  initial  condition  z_m+i  =  (e  ^  Hv  II  0O»  0,...,0). 

£ 

The  goal  of  the  following  analysis  is  to  show  that 


(3.33) 


-i( 0-£)h 

sup  e 
i<0 


Hz, 


II 


<  const  Hv  It 

t 

l 


which  implies 


(3.34) 


-i( a-e)h 

sup  e 
i<0 


|ui' 


<  const. 


-i( a- e)h 

sup  e 
i  <0 


Summing  up  the  geometric  series  in  (3.32),  we  obtain  the  estimate 


e_1<  0~e)hllzi  »  <  const. 


( i+m-1 ) (  o-  e) h 


(3.35) 


(<e 


(  0-e)h 


-  1)  II  [l 


+  1 


I+A(h, y) 


l 

■1  1 


(o-E)h  '  j 


We  thus  have  to  prove  estimates  of  the  following  form 


(3.36) 


fI+A(h, y) i 
^  (  o-  e  )  h  ' 


( A 
0 

u 


<  const . ,  k  6  N 


and 
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(3.37) 


( a-e)h 
e 


«(I 


I+Mh.y)  i 

(c-e)h  J 

e 


<  const. 


Both  will  follow  from  an  analysis  of  the  Jordan  form  of  A(h,y).  For  h,y 
small  enough  we  can  write 

(3.38)  A(  h , y) 

where 


(3.39) 

A( 0 , 0 ) 


This  matrix  has  the  characteristic  polynomial 

N  N  X  K 

(3.40)  p(p)  =  C(0,0)  l  K  -  p  -  C(0,0)  l  . 

*=1  *  *=1  VP 

N  -1 

Recalling  that  C(0,0)  =  (  £  — )  ,  it  is  an  easy  exercise  to  show  that  the 

H=1  AJl 

root  0  is  two  fold.  An  analysis  similar  to  that  given  for  (2.3)  shows  that 

all  remaining  zeros  are  real  and  negative. 

A  similar  calculation  shows  that  zero  is  also  a  double  eigenvalue  of 

A(0,0)  +  O(h)  (as  of  (3.38)).  Therefore  we  get  for  the  eigenvalues  of 
dh 

A(  h ,  y) 

(3.41)  (a)  P.j(h,Y)  =  ho(  1  ) ,  p  (h,y)  =  ho(  1 )  as  y  +  0 

(3.41 )  (b)  p  (h, y)  =  h( pi  +  o( 1 ) )  as  h,y  +  0,  i  =  3(1)(N+1) 


C(  0 , 0 )  l  -CfO.Oiy..  -0(0,0)^ 


-X. 


\K» 


=  h(—  A(  0 , 0  ) 
dh 


+  0(h)  +  0( y) ) 
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r 


where  <  0  for  i  =  3(1)(N+1)  such  that 

1  +  p^h,  y) 

{ o-e)h 
e 

holds.  Equality  in  (3.42)  only  holds  for  h  =  0.  Since  —  A(h,0)  is 

holomorphic  in  h  *  0  and  since  the  eigenvalue  cf  ^  A(h,0)  do  not  change 

d 

multiplicities  as  h  +  0  (the  negative  eigenvalues  of  — —  A(0,0)  are 

an 

distinct  and  0  is  a  double  eigenvalue  of  ry  A(h,0)).  There  is  a  matrix 

h 

G(h)  such  that  G(h),  G-‘(h)  are  holomorphic  in  h  =  0  and  J(h)  defined 
by 

A(h,0)  =  G(h)J(h)G~\h) 

is  the  Jordan  form  of  A(h,0)  (see  [2]).  Therefore  (3.36),  (3.37)  hold  for 
y  =  0.  A  simple  perturbation  argument  assures  (3.36),  (3.37)  for  y 
sufficiently  small.  Thus  for  h  sufficiently  small  and  K  >  0  sufficiently 
large,  we  have  proved  that 


(3.42) 


(3.43) 


-t  (  a-  e )  _ 

sup  e  |u^|  <  const  llvll  ^ 

t.<-K  i°° 

i 


The  solution  u^  can  be  continued  over  the  finite  interval  [ — K , 0 ]  and  by  a 
standard  stability  analysis  (see  [1])  we  obtain  (3.34). 

We  now  have  to  treat  the  case  t  >  0.  For  this,  we  rewrite  (3.21)  as 


( 3.44) 


1-1 

u  =  h  y  a( ( i- j )h) a  ,(h,-m)u. 

1  j=-m+ 1  io  3 

i-1 

+  h  y  a((i-j)h)a,  ,(h,-m)u.  +  v. 


where  it  is  assumed  that  tj  >  K  is  suf f iciently  large.  After  some 
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calculation,  we  get  from  (3.15),  (3.22) 


(3.45) 


a.  (h,-m)  *  D(h )  +  0  (h,-m) 

i ,  j  i  j 


where 


(3.46)  D(h) 


1  -ati 

-= - n -  ;  16.  .  (h,-m)  |  =  0(e  J),  t  >  K 

N  ^  0  1  r  J  J 


h  l 


*=  i  V  , 

e  -1 


It  is  therefore  natural  to  study  the  equation 


i-1 


(3.47) 


u  =  D(h)  l  a ( ( i- j ) h ) u ,  +  v  , 
j=I  : 


where  is  plus  the  first  sum  in  (3.44),  and  interpret  (3.44)  as 

perturbation  of  (3.47).  As  before  we  put 


(3.48) 


ji  V  "WV- 

9i =  V  I  e  ui 

j=I  : 


which  leads  to  the  difference  equation 


gi+1 


(3.49) 


-X  h  N  -X  h 

g  =  K  D(h)he  *  £  qJ  +  (e  Z  -  1)gZ 

*  j=1  1  x 


-*tK 

+  h  Ke  v. 
a  i 


Here  the  relation 


(3.50) 


N 

u,  =  D(h)  J  g?  +  v. 
i  .  .  x  x 

3=1 


has  been  used. 
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r 


■ 


IN  In 

Putting  g^  =  (g  ,...,g^),  e(h)  :  =  O^e  ,...,K^e  ),  we  obtain  the 

following  matrix  form  of  (3.49) 

(3.51)  g  =  (I  +  B(h))qi  +  v  Jn  »e( h)  . 

This  has  the  solution 


(3.52) 


i-1 


l 

j=I 


(I  +  B(h))1-j_1 


v  ,e(  h) 
3 


When  dealing  with  the  case  t^  <  0,  we  used  a  'redundant'  system  of  equations 
rather  than  an  analogue  of  (3.49).  The  reason  for  this  was  that  it  is  easier 
to  compute  the  characteristic  polynomial  of  the  matrix  of  the  'redundant' 
system.  In  the  'redundant'  ( (N+1 ) -dimensional  instead  of  N-dimensional )  form 
(3.49)  reads 


(3.53) 


N  -X.h  __  N  -X  h 

u  -  u  =  D(h)h  l  K  e  u  +  D(h)  £  (e  -  1)g 

1=  1  *  1= 1 


+  v 


i+1 


v. 

l 


5  i+1 


K  h  e  u.  + 

l  l 


( e 


-X*h 


Dgf 


When  we  write  this  in  matrix  form 

(3.54)  z  =  (I  +  A(h))zi  +  dj, 

~  1  N  ~  ~ 

(where  z^  =  (ut'?i* •• • '<Ji> *  ^  =  <vi+1  “  Vi'  0,...,0)),  we  see  immediately 
that  A(h)  is  the  same  matrix  as  (3.31),  exceot  that  C(h,y)  is  replaced 
by  D(h).  The  characteristic  polynomial  is 


(3.55) 


[<p)  =  D(h)h  l  K 


-v 


-  p  +  D(h)h  l  K 


i=1 


1=1 


-X.h  -X.h 

e  1  (.  1  -  II 

,  '*‘h  ,, 

p-  ( e  -  1 ) 


It  is  easily  verified  that  p  =  0  is  a  double  root.  Moreover,  since 
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(3.56)  D(h)  =  C(0,0)  +■  0(h)  , 

the  other  roots  are  small  perturbations  of  those  of  p(  p)  as  given  by  (3.40) 
and  have  therefore  negative  real  parts.  When  passing  from  A(h)  to  the  N- 
dimensional  matrix  B(h),  the  eigenvalues  obviously  remain  the  same,  except 
that  0  as  an  eigenvalue  of  B(h)  has  multiplicity  one  rather  than  two. 

Hence  there  is  a  matrix  E(h)  such  that  E(h),E-1(h)  are  continuous  for  h  e 
[0,hg]  and  the  Jordan  form  J(h)  of  B(h) 


(3.57) 

has  the  block  form 


J(h)  =  E_1(h)B(h)E(h) 


/ : . \ 


(3.58) 


J  (h)  = 


J  (h) 


where  the  (N-1  )  x  (N-l)-matrix  J_(h)  has  only  eigenvalues  with  negative 

real  parts.  The  continuity  of  E(h),E-1(h)  holds  since  B(h)  is  analytic 

n 

in  h  and  since  the  eigenvalues  of  ~  B(h)  are  distinct  even  for  h  =  0 

n 

( see  [2] ) . 

If  we  put  g^  =  E(h)w^,  (3.52)  yields 

(3.59)  w.  =  h  Y  (I  +  j(h))1"^”1  v.E-1(h)e(h)  . 

l  .  i 

j=I  J 

In  the  first  component  this  reads  in  particular 


( 3.60  ) 


=  h 


i- 1 

v 

/. 

3=1 


(E-1(h)e(h))  1 


From  this  we  obtain  the  following  estimate  for  t  -*•  »>: 
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(3.61  ) 


■i  •  i  1  i  ,  i(  o-e)h ,  ~  , 

lim  |w  |  <  const  sup  e  |v  | 

i  +«°  i  >1 


i(o-e)h.  1  1.  i(o-e)h.~ 

sup  e  |w  -  lim  w. (  <  const,  o ( 1 )  sup  e  |v. 

i>I  1  i>i  1 


2  N 

For  the  components  (w^,...,w  ),  where  only  eigenvalues  with  negative  real 
parts  are  involved,  an  analogous  estimate  follows  from  the  same  arguments  that 

l 

have  been  used  in  the  case  t^  <  0,  and  we  even  have  lim  w.  =  0  for 

1  ->CD 

l  >  1. 

Let  us  now  introduce  the  spaces 

A={f=(f.).!  lim  e  |f.|  =  0} 

■h,1  i  i-i  i>t  i 


Bh,I  *  <f  *  1  U”  't  ’ 

1  +00 


:  f  exists,  lim  e^°  |f.-f  I  =  0} 

OO  1  oo 

i>I 


and  the  operator 


(3.62) 


_  , .  cr-  e  cr-c 

Vh)  •’  Vi  +Bh,I 


which  is  defined  as  the  solution  operator  corresponding  to  (3.47),  i.e.  the 

OO  -s/  00 

operator  mapping  (v  )  to  (u^)^_^.  w*len  we  put 


(3.63) 


(Gu).  =  h  V  a((i-j)h)  s.  .<h,-m)u.  , 

1  j-I  ^  3 


(3.44)  can  be  rewritten  in  the  form 


(3.64) 


u.  =  P  (h)(Gu  +  v)  ,  i  >  I 
l  I  i 


It  follows  from  (3.61)  that  P-j(h)  is  a  bounded  operator.  Moreover,  (3.46) 


implies  that  G  has  small  norm.  Therefore,  I  -  PjthiG  is  a  nonsingular 

operator  from  F'  into  E,  and  the  norm  of  the  inverse  is  bounded 

n  f  I  n ,  I 

uniformly  with  respect  to  h,  t  and  e.  Therefore, 
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(3.65)  sup  e £)h|u.  -  lim  u  |  +  |lim  u  |  <  c Rv S 

i>0  1  i-*-«  1  i+»  1  Y*3- e 

h 


and,  summarizing,  we  obtain 


(3.66) 


Hu#  <  const  flv # 

Yo-e  „o-e 

*h  \ 


where  the  constant  is  independent  of  h,  t  and  e  e  [0,e  ],  where 

-rn  o 

eQ  <  a.  This  concludes  the  stability  proof  for  y  =  0. 


We  briefly  sketch  the  stability  proof  for  y  >  0.  Equation  (3.19)  now 
takes  the  form 

•a 

i-1 

a  ( (  i- 1 )  h )  f  2  — r — —  +  it. 

3 


Ui"Ui-1 


(3.67) 


=  -Hi(h,-m)ui  +  h  l  a((i-j)h)(2  — — —  +  1  ]u  .  +  v . 


j=-m+1 


y  (t.) 


i  >  -m 


where 


t  .  i-1 

-m  , .  ,  2 , 


y2(t.  ) 


(3.68) 


H^h.-m)  =  3  aft^sjy  (tjds  +  3h  £  a((i-j)h)  — — - 

j=-m+1  '  ‘ 


y  (t.) 


-  af(t.)»ya  1  ( t .  ) 

l  l 


With  D(h)  as  in  (3.46), 


-ot. 

0(  e  X) 


(3.69) 


Vh'-m)  =  W 


as  t .  ♦  «> 
i 


at . 


0(h)  +  0(e  )  as  t.  -*•  -<*> 

i 


We  can  therefore  use  a  similar  perturbation  approach  as  before,  i.e.,  for 
t^  +  “,  (3.67)  is  regarded  as  a  perturbation  of  the  problem 


(3.70) 


u .  -u  .  , 

l  1-1 


i-1 


‘  mh7  ui  +  3h  .2  ,  a((i~i)h)“j  +  \ 

j=-m+1  J 


This  can  be  rewritten  as  follows 


-3  3- 


(3.71) 


1.-2 

u^  =  +  5(h#y)h  £  a((l-j)h)u.  + 

j=-m+1  ^ 


where 


T(h,u»  -  *  3h.0»»>  +  g^jy)  .  6(h,„)  -  3(^.5^)  . 


We  substitute 


(3.72) 


.  i-2  -X.(t .  -t.) 

g.  =  h  ft,  e  u. 

i- 1  l  .  L  ..  3 

3=-m+ 1 


N  -  X  „h  . 
~  v  l 

'  l e 


^  V  ■»>  -v  *W  ‘V 

We  set  =  (^  **£  ^  =  (SOxuHv^  - 

(3.70)  is  equivalent  to  the  system 

(3.73)  =  (X  +  F(h,u))wi_1  +  ^ 


,.,0).  Then 


with 


■y(h ,  u)  -  1 


<5(h ,  y ) 


(3.74) 
F(h,y)  = 


N  -2X  h 

h  l  V 

2=1  * 


h 


-X^ 


I  “XMh 

\  »  V 


-X  h  -X  h 
e  (e  -  1) 


-X,h 

e  -  1 


-v  ■ 

e  ( e 


-v 
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As  before,  it  can  be  shown  that  the  characteristic  polynomial  of  D,F(0,y) 

h 

has  the  same  roots  as  (2.3),  except  for  the  fact  that  0  is  a  double  rather 
than  a  simple  root.  Moreover,  0  is  an  exact  eigenvalue  of  F(h,y). 

A  proof  analogous  to  the  one  for  y  =  0  shows  the  stability  for  y 
fixed  and  sufficiently  small  h.  For  the  limit  y  +  0,  a  different  argument 
is  needed.  When  we  substitute  in  (3.73) 


(3.75) 


~  5(h,y) 

u.  +  — — - f— -  z.  =  p. 

l  y  ( h ,  y )  - 1  l  *1 


Zi  =  qi 


we  obtain  a  system  of  difference  equations  of  the  form  (3.73)  with  F(h,y) 
substituted  by  a  matrix  of  the  following  form 


(3.76) 


F(h,y)  = 


y(h,y)  -  1  0(h) 


h-  A( 0 , 0 )  +  0(h2) 


Moreover,  an  estimate  of  the  form  -1  <  y(h,y)  -  1  <  -o)min(C.,C  — )  holds 

1  2  y 

where  u>  >  0.  Thus  for  small  y,  | y— 1 |  >>  h.  It  is  easy  to  conclude  from 
this  that  there  is  a  coordinate  transformation  close  to  the  identity  which 
transforms  F  to  the  form 


(3.77) 


y(h, y)  -  1  +  0(h) 


h  —  A( 0 , 0 )  +  o(h) 
dh 


Stability  follows  from  the  above  estimate  for  y  -  1  and  an  analysis  of  the 

eigenvalues  of  — —  A(0,0)  given  by  (3.39).  For  t  +  -<*>  a  similar  argument 
dh  i 

holds,  but  D(h)  is  to  be  replaced  by  a  different  constant  D(h,t  ).  From 

— m 

these  considerations  we  see  that 


(3.78) 


<  const. 


J 
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with  a  constant  independent  of  h,  -m,  y  and  e  e  [0,e  1,  <  a. 

It  is  practically  important  to  assure  stability  not  just  for  h 
sufficiently  small,  but  also  for  arbitrary  h.  Recall  that  linearizing 


-\h  1 
e  -  1  with  respect  to  h  is  only  justified  if  h  <<  —  .  The  matrix 

F(h,y)  for  arbitrary  h  has  the  same  form  as  for  h  small,  if  the 

-V  ,  -2  X ,  h 

e  -  1  i 

following  substitutions  are  made:  X.  +  - - - ,  K.  +  K.e  , 

l  n  ii 


U  *  V  + 


3h 


If  f  has  compact  support,  this  is  sufficient  to  ensure 


D(h) 

stability.  If  the  support  of  f  is  not  compact,  stability  for  arbitrary  h 


can  be  assured  if  the  following  modification  is  made:  In  (3.3)  the  integral 


t_  -m 

f  a(t.-s)ds  is  replaced  by  h  •  T  a(t,-t.).  With  this  modification 

J  -oo  i  r  ,L  1  1 

j  =  -  00 

the  term  0(h)  in  (3.69)  vanishes,  and  thus  the  matrix  for  the  linearized 
problem  is  asymptotically  equal  to  F(h,u)  both  for  -*•  <*•  and  t  * 

In  order  to  apply  Keller's  [3]  nonlinear  stability  concept,  it  is  further 
necessary  to  show  that  the  Frechet  derivatives  DgF^  m  are  uniformly 
Lipschitz  continuous  in  a  sphere 


S  =  {z  e  X°“E  |  II z  -  (y(t.)  -  1)”_  II  <  K }  . 

K  n  i  i  °° 

This  follows  from  a  fairly  trivial  calculation,  which  we  do  not  present  her.'. 

CO 

Using  the  fact  that  the  global  error  (y(t.)  -  y  )  is  estimated  bv  i 

11  1— 

constant  times  the  bound  for  the  local  error  (3.12),  we  obtain  the  following 
theorem: 

Theorem  3, 1 

The  discretization  scheme  (3.3),  (3.4)  has  a  unique  solution  for  a )  l 

J ,  1  A 

f  €  X  ,  where  a  is  as  of  Theorem  2.3,  this  solution  y  -  (v.l  <-in 

1  l  =  -  ■*> 

he  calculated  by  the  Newton  procedure  which  in  second  order  converuenf  from  i 


-  U  - 


sphere  of  starting  values  which  does  not  shrink  to  ?  as  h  +  0,  t  -*•  -°°, 

-m 

i  -►  0  and  the  convergence  estimate 

-e(t_  ) 

(3.79)  II  (y  -  y(  t ,  )  )  .  II  <  const  (h  +  o(  e  m  )  ) 

i  l  i=-«>  o-£ 

\ 

holds  for  h  sufficiently  small  and  |t_m|  sufficiently  large.  The  constant 

is  independent  of  h,  t_m,  y  e  [0,“],  £  e  [0,e],  <  o. 

This  implies  that  the  Newton  procedure  for  the  solution  of  (3.6)  can  be 

00 

safely  applied,  that  the  (y^K_  ^  do  not  exhibit  boundary  layer-like 
behavior  and  that 


(o-e)t.  -e|t_  | 

(3.80)  |y^  -  y(t^)|  <  const,  e  1(h  +  o(e  m  )),  t^  <  0 


~  e 1 t_m I 

(3.81)  |  lim  y.  -  y(«>)  |  <  const,  (h  +  o(e 

i+co 


) ) 


(3.82)  I  ( y .  -  lim  y.)  -  (y(t.)  -  y(<”)|  <  const,  e 

i+« 


-  (  a-  e )  t .  -  e  1 1  1 

1 ( h  +  o( e  “m  ) ) ,  t  >  0 


and  the  order  of  convergence  is  independent  of  ye 

Obviously  if  f  is  only  supported  on  [T.,T7],  then  the  term 

-elt_ml 

0(e  )  disappears  from  the  error  estimate  if  t  <  T^. 

The  discretization  we  used  was  derived  from  the  integral  equation.  In 
chapter  2  we  transformed  to  a  system  of  ordinary  differential  equations.  In 
fact,  up  to  terms  of  order  0(h),  our  discretization  method  corresponds  ‘..5  a 
discretization  scheme  for  the  ODE  system  (2.1).  Namely,  if  we  put 


’e,  i 


=  h 


j=-m+ 1 


V 


■X  Jt.-t .  ) 
l  i  ] 


Y3 


a 


l ,  -m 


0 
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0 


i  -A  (t  -t.) 

Vi  = h  .  I  V  3  y^  '  h 

3=-m+1 


3 


£,-m 


our  discretized  equation  reads  as  follows: 


-u( 


yi“yi-1 


N  -A.h 

)  =  l  e  (9.  t.Ji  -  h, 


£=1 


£,  i- 1  i  i, i-1 


(a) 


-  f.  yf  +  f  a(t. -s)ds(y  -  1) 
i  1  •'-»  l  i 


(3.83) 


q  -q  K  "Vh 

g£,i  g£,i-1  t  (e  -  1) 

(b)  h  ■  “ 2  +  h  Vi-1 

yi 

.  ,  -Ah 

.  .  Vi  '  h£, i- 1  ..  (e  1  -  D  , 

(c)  - - -  «  Vi  +  - - -  Vi-1 


By  calculating  g  h  .  for  £  =  1(1)N  from  (3.83)(b),  (3.83 )(c)  and  by 
X.  /  1  X-  #  X 

inserting  these  quantities  into  (3. 83) (a)  an  equation  of  the  fon>  (3.6)  is 
obtained  (in  each  time  step).  Theorem  3.1  now  implies  that  the  root  of  this 
equation  can  be  safely  obtained  by  the  Newton  procedure  which  is  second  order 
accurate  from  a  sphere  of  starting  values  whose  radius  is  independent  of  h, 
t  ,  u  e  [0,”)  and  i  >  -m. 

This  provides  us  with  a  very  efficient  method  to  solve  the  approximating 
problems  and  Theorem  3.1  makes  sure  that  the  qualitative  properties  of  the 
solution  of  (1.1)  carry  over  to  the  approximate  solutions. 

The  exponential  decay  of  the  solution  encourages  one  to  attempt  using 
variable  mesh  sizes  of  the  form 

a  -alt  I 

(3.84)  h  =  he 
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A 

It  can  be  expected  that  convergence  of  the  order  one  in  h,  i.e.  the 

A  _o|t 

A  .  —  Jfl  *  00 

estimate  would  be  0(h)  +  o[e  J  ,  would  follow  in  i.  but  the 

exponential  decay  property  of  the  approximate  solution  would  be  lost.  In  the 
case  of  boundary  value  problems  for  ordinary  differential  equations  on 
infinite  intervals  this  has  been  shown  in  [7] . 

A  furthei  problem  that  should  be  mentioned  is  which  higher  order 
discretization  schemes  could  be  employed.  It  is  fairly  clear  from  our 
analysis  that  polynomial  collocation  methods  using  Radau  points  (see  [12]) 
could  be  used  and  the  exponential  decay  property  and  the  uniform  convergence 
as  u  -*•  0  would  be  recovered. 
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4 


Numerical  Results 


For  the  computations  we  used  the  kernel  a(u) 


followinq  constants  K.  and  X 

i  i 


a  -x.u 

T  K.e  1  with  the 
L .  i 


i 

-1 

X . ( sec  ) 
i 

Ki(Nm~2sec_1 ) 

i 

10-3 

— fc 

X 

_* 

o 

1 

u> 

2 

10~2 

o 

o 

X 

CO 

• 

r— 

3 

10'1 

2 

1.89  x  1(T 

4 

1 

9.8  x  103 

5 

10 

2.67  x  105 

6 

102 

5.86  x  106 

7 

103 

9.48  x  107 

8  1 

104 

1.29  x  109 

These  numbers  were  obtained  by  Laun  [4]  from  an  experimental  fit  for  a 

polyethylene  melt  at  150°C,  which  he  calls  "Melt  1". 

The  parameter  u  is  physically  identified  as  three  times  the  Newtonian 

contribution  to  the  viscosity.  Experimental  values  are  not  available,  and 

theoretically  u  is  either  a  solvent  viscosity  (for  polymer  solutions)  or  it 

results  from  fractions  of  low  molecular  weight  (for  melts).  The  value  of  u 

has  to  be  compared  to  the  viscosity  resulting  from  the  memory,  which,  for 

ft  _2 

constant  shear  rate,  is  given  by  \  K.  X.  *  50000  Nm'^sec.  would  expe<-> 

‘  li 

i=  1 

U  to  influence  the  solution  significantly  only  if  it  exceeds  this  value. 
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This  is  verified  by  our  computations.  In  the  plots  the  scale  for  y  is  on 

the  left,  the  scale  for  f  is  on  the  right,  y  is  measured  in  multiples  of 

the  length  (for  the  filament,  a  =  2),  or,  respectively,  the  thickness  (for 

the  sheet,  a  =  j)  at  t  =  f  denotes  the  force  acting  on  the  ends  of  the 

filament  or  the  edges  of  the  sheet  divided  by  the  cross-sectional  area  in  the 

N 

undefonned  state  at  t  =  -°°,  f  is  expressed  in  — .  The  time  is  measured  in 

m 


seconds,  f  is  always  plotted  by  dashed  lines,  y  by  full  lines. 

All  plots  except  figures  7,  8  were  made  for  a  =  2,  the  case  of  the 
filament.  In  figures  1-11  (except  6),  the  force  f  is  of  the  form 

f  0  Itl  >  a2 


f(t)  = 


max 


exp  {a 


1 


2  2 
V1 


with  a. 


1 

-  “a  -  °* 
a2 


Such  an  f  is  in  C  (*, R)  and  has  the  compact  support 


[-a2,a2] . 

The  parameter  y  is  zero  in  figures  1-9.  In  figures  1  -  5  we  have 
chosen  various  values  of  f  max,  a^  and  a^,  as  can  be  seen  from  the 
diagrams.  The  calculations  were  done  for  larger  time  intervals  than  the 
plots,  thus  yielding  approximations  for  y(»).  For  figures  1-5,  the 
approximate  values  of  y(  <*>)  are  as  follows: 


fig. 

1 

2 

3 

4 

* 

5 

y  ( °°) 

1.07 

1.15 

1.11 

1. 25 

1.13 

* ( in  this  figure  supp  f  is  different  from  the  previous  ones) 
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Figure  4 


Figure  6 


r 


These  numbers  indicate  that,  roughly  speaking,  the  value  of  y(«)  -  1  is 

rOO 

proportional  to  I  f(t)dt.  This  would  in  fact  be  exact  for  the  linearized 

i  —00 

equation. 

In  figure  6  an  oscillating  force  was  chosen.  It  is  observed  that  the 
solution  y  "follows"  the  oscillations  with  a  certain  delay. 

The  figures  7  and  8  illustrate  the  case  of  the  sheet  (a  =  “■)  •  Here 
-f  is  plotted  rather  than  f.  The  results  are  qualitatively  similar  to  those 
in  figures  1-5,  but  now  we  have  y  <  1  instead  of  y  >  1. 

In  figures  9-15,  we  have  again  ex  =  2.  In  figures  9  -  11,  we  have 
chosen  the  same  f  (f  max  =  40000,  Sq  =  1,  a1  =  aj  =20)  and  computed 


solutions  for  different  values  of 

P 

p* 

y  1 

y  max 

y  ( °°) 

0 

1 

i 

103 

3.3 

105 

89 

3.4 

106 

7.7  I 

4.5 

For  u  <  10000,  no  significant  change  was  observed.  For  larger  u,  the 
effect  on  the  maximal  elongation  seems  to  be  more  pronounced  than  the  effect 
on  the  final  length.  Recalling  the  fact  that  \i  =  1h000  would  correspond  to 
a  viscosity  3  x  10^  as  large  as  that  of  water,  it  seems  conceivable  that  for 
fluids  like  "Melt  1"  u  can  be  neglected. 

The  numbers  for  y(®)  are  interesting  in  comparison  with  the  results  of 
Lodqe ,  McLeod  and  Nohel  [6).  They  showed  that  y(°°)  increases  with  p,  if 

the  history  of  y  for  t  <  0  is  kept  fixed.  Our  numbers  show  the  same 
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tendency,  eventually,  however,  y(°°)  has  to  decrease,  since  for  y  =  *>  we 


have  y  =  const,  and  thus  y(“)  =  1.  We  see  from  this  that,  for  fixed  f, 
y(®)  is  not  a  monotone  function  of  y. 

For  figures  12  -  15,  a  discontinuous  force  given  by 


f(t) 


JO  t  >  0 

40000  exp(^)  t  <  0 


was  used.  Since  in  this  case  the  filament  recovers  freely  for  t  >  0,  we  are 
studying  the  same  situation  as  lodge,  McLeod  and  Nohel  [61,  but  we  prescribe 
the  force  rather  than  the  history  of  y  for  t  <  0.  By  considering  the 
intervals  t  <  0  and  t  >  0  separately,  we  can  easily  modify  the  existence 
and  convergence  theory  of  the  previous  chapters  for  the  present  case . 

CO 

However,  the  solution  does  not  depend  continuously  on  y  in  the  L  -norm  as 
y  ♦  0.  This  is  because  for  y  =  0  the  solution  is  discontinuous  at  t  =  0. 
The  following  table  illustrates  the  dependence  of  the  maximal  elongation  on  y 


u 

y  max 

0 

140 

2  x  10S 

58 

3.5  x  105 

16 

10* 

1.6 
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ABSTRACT  (continued) 


solution  which  approaches  a  given  limit  as  t  -*■  moreover,  the  r  olution 

also  has  a  limit  as  t- -*■  +«.  A  numerical  scheme  is  analyzed  and  convergence 
uniformly  in  t  is  established.  Particular  attention  is  paid  to  the 
dependence  of  solutions  on  a  parameter  y,  which  corresponds  to  a  Newtonian 
contribution  to  the  viscosity.  It  is  proved  that  solutions  converge  uniformly 
in  t  as  y  -*  0,  and  that  the  convergence  of  the  numerical  scheme  is  also 
uniform  in  y. 
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